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ABSTRACT 

Context. Galactic open clusters have been long recognized as one of the best tools to investigate the chemical content of Galactic disk 
and its time evolution. In the last decade, many efforts have been directed to chemically characterize the old and intermediate age 
population; surprisingly, the chemical content of the younger and close counterpart remains largely undetermined. 
Aims. In this paper we present the abundance analysis of a sample of 15 G/K members of the young pre-main sequence clusters 
IC 2602 and IC 2391. Along with IC 4665, these are the first pre-main sequence clusters for which a detailed abundance determination 
has been carried out so far. 

Methods. We analyzed high-resolution, high S/N spectra acquired with different instruments (UVES and CASPEC at ESO, and the 
echelle spectrograph at CTIO), using MOOG and equivalent width measurements. Along with metallicity ([Fe/H]), we measured 
Na i, Si i, Ca i, Ti i and Ti n, and Ni i abundances. Stars cooler than ~ 5500 show lower Ca i, Ti i, and Na i than warmer stars. By 
determining Ti n abundances, we show that, at least for Ti, this effect is due to NLTE and over-ionization. 

Results. We find average metallicities [Fe/H] =0±0.01 and [Fe/H]=-0.01 + 0.02 for IC 2602 and IC 2391, respectively. All the 
[X/Fe] ratios show a solar composition; the accurate measurements allow us to exclude the presence of star-to-star scatter among the 
members. 
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1. Introduction 

Open clusters (OCs) are widely recognized as excellent tracers 
of the chemical properties of the Galactic disk and of its evolu- 
tion with time. In recent years, several studies have been devoted 
to the determination of the chemical composition of the older 
and intermediate age population, i.e. OCs with ages greater than 
~ 500 Myr, with the main goal of deriving the radial metallicity 
gradient across the Galactic disk (e.g., Carraro et al. 2004; Yong 
et al. 2005; Friel 2006; Randich et al. |2006l Carretta et al. |2007l 
Sestito et al. 2007, 2008; Pace et al. 120081 . On the other hand, 
less attention has been paid to the abundances of younger clus- 
ters and, in particular, the so-called pre-main sequence (PMS) 
clusters (age ~10-50 Myr). As we will discuss below, knowledge 
of the abundance pattern in these clusters is instead important to 
address different astrophysical issues, such as the evolution of 
debris disks, the common origin of young open clusters and as- 
sociations, and the chemical evolution of the solar vicinity. 

First, the age interval defined by PMS clusters is critical to 
trace the early stages of evolution of dusty debris disks around 
main-sequence stars. The fraction of debris disks in clusters and 
the level of 24 fim excess emission rise for stars with ages be- 
tween 5-10 Myr, reach a peak around 20-30 Myr, and then start 
decaying inversely with age (e.g., Currie et al. 2008). Planet for- 
mation models suggest that terrestrial planets also reach their 
final mass by 10-30 Myr; hence, a good characterization of de- 
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bris disks in this critical age range provides an indirect tool to 
study planet formation and evolution (e.g Siegler et al. 2007). 
Several studies have focused on the detection and characteriza- 
tion of debris disks both in the field and in open clusters and 
on the existence of possible correlations between the frequency 
of dusty disks and the fundamental properties of the parent- 
star (spectral type, mass, and luminosity for instance). Given 
the metal-rich nature of stars with giant planets (e.g. Santos et 
al. 2004), the correlation between metallicity and frequency of 
debris disks has been (partially) investigated: the lack of such 
correlation (Beichman et al. 2006; Bryden et al. 2006; Trilling 
et al. 2008) might reflect the different formation histories of gi- 
ant planets and debris disks. On the other hand, the latter stud- 
ies have focused on old solar-type stars in the field (ages ~ few 
Gyr) and the younger population has not been involved at all: 
there is instead the possibility of an initial correlation between 
dust production and metallicity around young stars, but this re- 
lationship might have disappeared as the stars age (Bryden et al. 
2006; Dominik & Decin 2003). PMS clusters and their young 
population represent a unique tracer to address this topic. In par- 
ticular, the determination of the metallicity of a large sample of 
PMS clusters with available information on debris disks would 
allow investigation of whether a correlation between metallicity 
and frequency of debris disks is indeed present at young ages. 

In a different context, the identification of small close-by as- 
sociations (within ~ 100 pc from the Sun) along with the knowl- 
edge of several complexes (e.g. the Sco-Cen association and the 
Tucana/Horologium) a new impulse has been given to the in- 
vestigation of the so-called superclusters. The concept of such 
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stellar aggregates with parallel space motions and located in the 
solar surroundings was introduced by Eggen more than thirty 
years ago; one of these superclusters is the Local Association 
(Eggen 1975, 1983a,b): with the Pleiades, other clusters like a 
Persei, 5 Lyrae, NGC 1039, NGC 2516, and the Sco-Cen OB 
association, also IC 2602 would be contained within it (Ortega 
et al. 2007). The secure knowledge of the chemical composi- 
tion and the detection of a common/different abundance pattern 
is critical to confirm or not a common origin, independently of 
dynamical/kinematical information. 

Finally, young clusters in the solar vicinity allow us to chem- 
ically characterize our immediate surroundings in the disk, the 
solar neighborhood, at the present time; at variance with old field 
stars that may have moved from their original birthplace, these 
clusters are very likely born close to where they are now and thus 
their abundance pattern should be representative of the present- 
day composition of the solar vicinity. Therefore, the comparison 
of the abundance pattern with that of older stars, including the 
Sun, allows us to place constraints on models of Galactic evolu- 
tion in the solar neighborhood. 

To our knowledge, only one PMS cluster has an available 
metallicity determination. Shen et al. (2005) presented an abun- 
dance analysis of 18 F-K dwarf members of the ~30 Myr old 
open cluster IC 4665. They obtained a very close to solar metal- 
licity ([Fe/H]=-0.03) with no dispersion among the members 
(standard deviation 0.04 dex), within the observational uncer- 
tainties. 

With this background in mind, we started a project aimed to 
derive the chemical composition in young clusters. As part of 
this effort, we present here an abundance analysis of G/K-type 
members of the two PMS clusters IC 2602 and IC 239 1 . 

The paper is organized as follows: in Sect. [2] we present the 
target clusters, observations, and the data reduction; the analysis 
is described in Sect. [3] while the results are presented in Sect. |4j 
the scientific implications are discussed in Sect. [5] followed by 
the conclusions (Sect. |6j. 

2. Target clusters, observations, and data reduction 

IC 2602 and IC 2391 are both PMS clusters with estimated ages 
of ~ 30 Myr (Stauffer et al. 1997) and ~ 55 Myr (Barrado y 
Navascues et al. 1999, 2004), respectively. The age of IC 2602 
has been estimated by classical isochrone fitting, while the age of 
IC 2391 has been derived using the lithium depletion boundary 
method. 

Due to both their age and distance from the Sun, several 
studies have focused on these two clusters, aiming to deter- 
mine different fundamental characteristics, such as membership, 
the lithium distribution among G/K-type stars, the X-ray prop- 
erties, rotational velocities, chromospheric activity, and initial 
mass function (Randich et al. 1995, 1997, 2001; Stauffer et al. 
19 89, 19 97; Prosser et al. 1996; Patten & Simo n 1996; Spezzi et 
al. I2009I I. More recently, Siegler et al. ( 120071 1 presented 24 pm 
Spitzer MIPS observations of the central field of IC 2391, which 
allowed them to to determine the fraction of debris disks around 
stars across a wide range of spectral types, from A to K. The 
comparison with field stars and other clusters indicates that the 
fraction of disks among A-type stars is lower than expected. 

Surprisingly, no extensive abundance study has been carried 
out for IC 2602 and IC 2391; Randich et al. (2001) derived the 
iron content ([Fe/H]) for a small sub-sample of stars in both clus- 
ters, finding close-to-solar metallicities with rather large stan- 
dard deviation around the average. Our goal here is twofold: first, 
based on higher quality spectra, we will derive a more accurate 



[Fe/H] for the two clusters. Second, besides iron, we will also 
determine abundances of other elements. Our sample includes 
eight and seven members of IC 2602 and IC 2391, respectively, 
for a total of 15 stars. From the initial and larger sample of con- 
firmed members in Randich et al. (2001), we selected only high- 
quality spectra of G- and K-type stars with rotational velocities 
(vsin i) below ~ 20 km/sec, since the analysis is much less reli- 
able for stars with higher v sin i values. 

The observations were carried out at the European Southern 
Observatory (ESO) and at the Cerro-Tololo Inter-American 
Observatory (CTIO). New observations were obtained for 12 
stars using UVES on VLT/UT2 (Dekker et al. 2000) during three 
nights of observations on 16, 17 and 18 February of 2001; a slit 
of 0.8", providing a resolving power of R~40000, and the CD4 
cross-disperser (5770A - 9420 A) were used. Only the blue part 
of the spectrum (5750 - 7450 A) was employed for the analysis, 
for consistency with the analysis of the spectra obtained with 
other spectrographs (see below). For the remaining three stars 
we used the spectra already published by Stauffer et al. (1997) 
and Randich et al. (1997 2001). Specifically, two stars were 
observed at CTIO with the 4m telescope in conjunction with 
the Red Long-Camera echelle spectrograph (RLCeS) and a 3 1 .6 
lines mm -1 . A 120 pm slit width (0.8 arcsec on the sky) and a 
Tektronix 2048x2048 CCD provided a resolution of R -43800 
and a spectral coverage between 5800 and 8200 A. Finally, we 
also added one star, IC2602-R95, whose spectrum was acquired 
on April 1994 with the CASPEC echelle spectrograph mounted 
at the 3.6m ESO telescope. A 31.6 lines mrrT 1 grating together 
with the red cross disperser was used. The slit width was 300 /mi 
(2.25 arcsec on the sky), resulting in a resolving power of R~18 
000 and a spectral coverage from 5585 A to 8400 A. 

Information on the target stars is provided in Table [1] while 
their spectra are shown in Fig. Q] where we plot a 60 A wide 
range from 6700 A to 6760 A. For details on the reduction of 
CTIO and CASPEC spectra we refer to Stauffer et al. (1997) and 
Randich et al. (2001), respectively. The UVES data were reduced 
using the UVES pipeline Data Reduction Software (Modigliani 
et al. 2004). The software is designed to automatically reduce 
point source object spectra from raw frames to order extraction 
and final order merging, by following the standard procedure: 
background subtraction, flat-field correction, order extraction, 
sky subtraction and wavelength calibration. 



3. Analysis 

3.1. Line list and equivalent widths 

The abundance analysis was carried out using the driver abfind 
of the spectral code MOOG (Sneden, 1973 -2000 version) and a 
grid of 1-D model atmospheres by Kurucz (1993). Radiative and 
Stark broadening are treated in a standard way; for the collisional 
damping we used the classical Unsold (1955) approximation. As 
discussed by Paulson et al. (2003) this choice should not greatly 
affect the differential analysis with respect to the Sun. LTE abun- 
dance values were derived by means of the equivalent widths 
(EWs) which were measured using a Gaussian fitting procedure 
with the IRAF task SPLOT . 

For the analysis we adopted as a basis the line list by Randich 
et al. (2006), but i) we selected only lines at A greater than 
5750 A; ii) we integrated the list with a few lines from Sestito et 
al. (f2006), choosing the lines included in our wavelength range 
and for which Sestito et al. obtained solar abundances in agree- 
ment with the standard values. We refer to those two papers for 
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Table 1. Sample stars: the first eight stars (R1-W79) are members of IC 2602, while the remaining seven belong to IC 2391. 
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star 


inst. 


exp.time 


S/N 


V 


B-V 


T c ff 


v sin i 






(s) 








(K) 


(km/s) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Rl 


TTVFS 

U V I -. i-) 


J AO / \J 




1 1 fS 

1 1 .VJ 




JUJU 


i 1U 


R14 


UVES 


5x870 


230 


11.6 


0.87 


5150 


13 


R15 


UVES 


10x870 


200 


11.7 


0.93 


4810 


7 


R66 


UVES 


4x860 


220 


11.1 


0.68 


5560 


12 


R70 


UVES 


5x900 


250 


10.9 


0.69 


5700 


12 


R92 


UVES 


7x300 


200 


10.3 


0.69 


5630 


14 


R95 


CASPEC 


1200 


150 


11.7 


0.87 


5020 


12 


W79 


RLCeS 


900 


100 


11.6 


0.83 


5260 


8 


VXR3 


UVES 


5x900 


260 


10.9 




5590 


10 


VXR31 


UVES 


7x900 


280 


11.2 




5630 


17 


VXR67 


UVES 


8x870 


130 


11.7 




4750 


8 


VXR70 


UVES 


5x900 


260 


10.8 


0.64 


5557 


17 


VXR72 


UVES 


8x900 


250 


11.5 


0.73 


5260 


15 


VXR76a 


RLCeS 


900 


70 


12.8 


1.05 


4400 


8 


SHJM2 


UVES 


7x300 


270 


10.3 


0.57 


5970 


<15 



(1) Identifier from Randich et al. ( 1995 1 for IC 2602 stars with the exception of W79 (Whiteoak 1961), and from Patten and Simon ( 1996 1 for 
IC 2391 members; (2) instruments (UVES, CASPEC and Red Long-Camera echelle spectrograph -RLCeS); (3) exposure time and number of 
exposures; (4) S/N r atios a t 6707.78 A; (5) V magnitude from Randich et al. (2001); (6) B-V colors from Randich et al.( 2001); (7) initial T eff 
from Randich et al. ( [20011 , and (8) vsin (' from Stauffer et al. (1997). 



Table 2. Solar abundances derived from our analysis for five 
elements; in the last column we report the standard values by 
Anders & Grevesse ( 119891 ). 



Element log n(X) log n(X) A G89 



Nai 


6.33+0.02 


6.33 


Si i 


7.56±0.01 


7.55 


Cai 


6.35±0.02 


6.36 


Tii 


4.97±0.01 


4.99 


Tin 


4.92±0.10 


4.99 


Nil 


6.26±0.02 


6.25 



further details on atomic parameters (log gf) and their sources; 
our line list for all the elements with EWs measured both in the 
solar spectrum (see next section) and in the spectra of our sample 
stars is available as electronic table. 

3.2. Solar analysis 

Our analysis was carried out differentially with respect to the 
Sun; as a consequence, most of the systematic errors due to e.g. 
atomic parameters should be canceled out. We analyzed the solar 
spectrum obtained with UVES (see Randich et al. 2006) using 
the same code, line list, and model atmospheres employed for 
our sample stars. The following stellar parameters were adopted 
for the Sun: T eff =5770 K, logg=4.44, £=1.1, and [M/H]=0. We 
obtained log n(Fe I) = 7.52+0.02 and log n(Fe II)=7.52 ±0.03. 
The abundance values derived for the other elements are listed 
in Table [2] together with the the values by Anders & Grevesse 
(1989) that were used as input in MOOG. As the table shows, 
the agreement is very good for all elements. 

3.3. Stellar parameters 

Initial effective temperatures (T e ff) were retrieved from Randich 
et al. ( 120011) who, in turn, derived them as using both B-V and 
V-I colors and the calibrations of Soderblom et al. ( 19933 and 



Bessell ( 119911 ), respectively -see Randich et al. (120011) for de- 
tails. 

We formally derived final temperatures by removing the 
trends between log n(Fe) and the excitation potential (x), after 
applying a 1-cr clipping criterion to the initial line list. In most 
cases, however, initial temperatures were already close to the fi- 
nal ones and only small adjustments were needed. This is shown 
in Fig. |2]where we compare our final T e ff values with those de- 
rived by Randich et al. (2001): as one can see, the agreement is 
good. 

For surface gravities, we adopted as initial values logg=4.5 
for all our sample stars. For the warmest stars in the sam- 
ple (T e ff > 5550 K) and for the very slow rotators (vsin/ < 
15 km/sec) we were able to check this assumption by impos- 
ing that the difference between the average log n(Fe I) and 
log n(Fe II) should be below 0.05 dex. More in detail, we 
were able to measure the Fe II for the stars R66, R70 and R92 
(IC 2602) and VXR3 and SHJM2 (IC 2391): the final value that 
we derived for all these stars is logg=4.45 and very agrees well 
with the initial assumption. 

Initial microturbulence values (£) were estimated using the 
Nissen's relationship (1981): £ = 3.2 x 10~ 4 (T eff - 6390) - 
1.3(logg - 4.16) + 1.7 km/sec. Final values were retrieved by 
zeroing the slope between iron abundances and EWs in MOOG. 



3.4. Errors 

As usual, the final abundances are affected by random (internal) 
and systematic (external) uncertainties. 

Internal errors due to uncertainties in oscillator strengths 
(log gf) should cancel out, since our analysis is strictly differ- 
ential with respect to the Sun. Random errors due to EW mea- 
surements are instead well represented by the standard devia- 
tion from the mean abundance based on the whole set of lines 
(a"i in [Fe/H] -see Table [4]). The errors (o"i) in [X/Fe] listed in 
Table |4] were, in turn, obtained by quadratically adding the rms 
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6720 6740 
Wavelength {A} 



.2 3 




6730 6740 
Wavelength (A) 



Fig. 1. Spectra of all the sample stars in a wavelength window 
of 60 A. Members of IC 2602 and IC 2391 are shown in the 
upper and lower panels, respectively. Several spectral features 
are marked. 



for [Fe/H] and rms for [X/H]. 

Final abundances are also affected by the random errors due to 
uncertainties in the adopted set of stellar parameters, namely 
Teff, logg, and These sources of error were estimated by 
varying one parameter at time, while keeping the other ones un- 
changed; then we took the quadratic sum of these three contribu- 
tions in order to obtain cr 2 , for [Fe/H] and for all [X/Fe] ratios. 
We adopted random errors of +60 K and +0.2 km/sec in T e ff and 
respectively: variations in T e ff larger than 60 K would have 
introduced a trend in log n(Fe) and^-; similarly, variations in mi- 
croturbulence velocities larger than 0.2 kms ~' would result in 
spurious trends in abundances versus EWs. For the five stars for 
which we could optimize the surface gravities, the uncertainties 
in logg were estimated by varying this quantity until the differ- 
ence between log n(Fe i) and log n(Fe n) is larger than 0.05 dex, 
i.e. the ionization equilibrium condition is no longer verified; 
this is achieved for a variation of 0.15 dex (R66 and VXR3) and 
0.20 (R70, R92, SHJM2) dex. For all the other stars, we adopted 
a conservative value of 0.25 dex. In Table[3]we list the errors in 
[Fe/H] and [X/Fe] ratios due to uncertainties in stellar parame- 



5000 5500 

T cff {this study) 



Fig. 2. Comparison of the effective temperatures of our sample 
stars with those derived by Randich et al. (1200 11 1; the dashed line 
is the 1:1 relationship. 



ters for one of the coolest stars and for the warmest one in the 
sample. 

The total random errors (<x tot ) were then calculated by 
quadratically adding cr\ and cr 2 ; the standard deviation of the 
cluster mean, listed in Table [4] is always smaller than the cr tot 
for all our sample stars and all the elements, suggesting that our 
error estimates are conservative. 

A global estimate of external errors due to the model at- 
mospheres and to the code can be obtained by analyzing one 
star whose chemical content is well known in the literature; to 
this purpose, we derived abundances for the Hyades member 
vbl87 whose spectrum was acquired with the same instrument 
(UVES -see Randich et al. 2007). The metallicity that we derived 
([Fe/H]=0. 14+0.04) is in excellent agreement with all previous 
determinations (e.g. Paulson et al. 2003), also [Si/Fe], [Ca/Fe] 
and [Ti/Fe] well agree within the uncertainties (see Table|4]i. We 
conclude that our analysis should not be affected by major sys- 
tematic errors. 



4. Results 

Our results for [Fe/H] and all [X/Fe] ratios are reported in 
Table [4] along with errors due to EW measurements and stel- 
lar parameters, namely cr\ and <T2. These final abundances were 
retrieved by applying a 1-cr clipping to the initial line list for the 
iron and for the other elements with a not too small number of 
lines (such as Ca 1 and Ti 1). Adopted stellar parameters are also 
listed in Table g] 

As a verification of our method, we show in Fig.|3]the [X/Fe] 
ratios as a function of T e ff for all our sample stars: as one can 
see, for [Fe/H], [Si/Fe] and [Ni/Fe] no obvious correlation is 
present in the diagram; on the other hand, sodium, calcium and 
titanium abundances are significantly lower than average for the 
coolest stars in the sample (R15, VXR67, and VXR76a). Trends 
of [X/Fe] vs. T e ff for various elements were already found for old 



Table 4. Derived abundances. We list again the identifier of each star (1), final stellar parameters (Cols. 2, 3, 4), [Fe/H] values (5), and the [X/Fe] ratios (Cols. 6, 7, 8, 9, 10). The 
errors (cr\, cr 2 ) are also reported. Average abundances for iron and the other elements are also indicated. 



Star 




logg 




[Fe/H] 


[Na/Fe] 


[Si/Fe] 


[Ca/Fe] 


[Tii/Fe] 


[Tin/Fe] 


[Ni/Fe] 




(K) 




(km/sec) 
















P 1 

K 1 


JUjU 




1 1 ^ 




n 1 A-i-H nA-i-n ha 

— U. 14±U.U i +±U.UO 


n rn-uO n^-i-H no 
u.uz±u.uj±u.uy 


n m_i_n ha 
— u.uz±u.u / ±u.uo 


— u.uo±u.uj±u.uo 


n ooj_n n/i_i_n 1 1 


— U.Ul±U.UD±U.Uo 


R14 


5100 


4.5 


1.1 


0±0.05±0.05 


-0.07±0.06±0.06 


0.05±0.06±0.09 


0.01±0.06±0.06 


-0.02±0.07±0.06 


0.11±0.04±0.11 


-0.01±0.08±0.06 


R15 


4770 


4.5 


0.7 


0±0.06±0.06 


-0.35±0.07±0.08 


0.03±0.07±0.09 


-0.2±0.05±0.07 


-0.4±0.05±0.06 


0.30±0.05±0.13 


0.02±0.05±0.08 


R66 


5590 


4.45 


1.15 


-0.01±0.03±0.06 


-0.03±0.04±0.05 


0.02±0.05±0.05 


0±0.05±0.02 


-0.03±0.04±0.04 


0.01±0.08±0.09 


-0.03±0.05±0.04 


R70 


5760 


4.45 


1.1 


0±0.05±0.06 


0±0.06±0.05 


-0.04±0.07±0.06 


0.03±0.1±0.04 


0±0.05±0.03 


0±0.04±0.08 


-0.04±0.08±0.05 


R92 


5630 


4.45 


1.2 


-0.02±0.04±0.06 


0.02±0.04±0.05 


-0.02±0.06±0.05 


0.02±0.05±0.03 


0.01±0.07±0.02 


0.02±0.01±0.08 


-0.02±0.06±0.04 


R95 


5020 


4.5 


0.8 


0.02±0.03±0.05 


-0.13±0.05±0.06 


-0.02±0.06±0.07 


0.04±0.06±0.07 


-0.07±0.08±0.07 


0.24±0.04±0.11 


-0.05±0.09±0.07 


W79 


5260 


4.5 


1.1 


-0.02±0.08±0.06 


-0.12±0.09±0.06 


0.01±0.08±0.07 


0±0.08±0.07 


0.02±0.08±0.05 


0.12±0.01±0.10 


-0.03±0.12±0.04 


IC 2602 








0±0.01 


-0.10+0.12 


0+0.03 


-0.02±0.08 


-0.06±0.14 


0.13±0.12 


-0.02±0.02 












-0.07 ± 0.06" 




0.01 ± 0.02" 


-0.01 ±0.04° 
















0+0.02 c 












VXR3 


5590 


4.45 


1.15 


0+0.04+0.07 


-0.02±0.05±0.05 


0.02±0.06±0.05 


0.01±0.05±0.02 


0.00±0.06±0.03 


0.03±0.03±0.09 


-0.02±0.06+0.04 


VXR31 


5630 


4.5 


1.2 


-0.02±0.04±0.07 


0.02±0.05±0.05 


0.02±0.05±0.05 


0.01±0.04±0.01 


0.03±0.08±0.02 


0.07±0.03±0.09 


0.02±0.06±0.02 


VXR67 


4680 


4.5 


0.7 


0.01±0.05±0.06 


-0.31+0.09±0.07 


0+0.09+0.09 


-0.21±0.11±0.08 


-0.38±0.04±0.06 


0.25±0.07±0.13 


-0.02±0.06±0.06 


VXR70 


5557 


4.5 


1.15 


-0.03±0.05±0.07 


-0.06±0.07±0.05 


0.04±0.06±0.07 


0.03±0.06±0.03 


-0.03±0.07±0.04 


0.09±0.04±0.10 


0.02±0.05±0.05 


VXR72 


5257 


4.5 


1.1 


-0.02±0.03±0.06 


-0.12±0.05±0.06 


0.01±0.04±0.08 


0.03±0.08±0.05 


0.02±0.05±0.05 


0.17±0.05±0.11 


0.01±0.04±0.05 


VXR76a 


4343 


4.5 


1.2 


0±0.06±0.07 


-0.34±0.08±0.07 


0.06±0.09±0.08 


-0.38±0.12±0.08 


-0.35±0.16±0.06 


0.37±0.06±0.13 


0.02±0.09±0.08 


SHJM2 


5970 


4.45 


1.2 


0+0.03+0.05 


-0.03±0.04±0.04 


-0.01±0.06±0.06 


0.02±0.04±0.04 


0.01±0.04±0.01 


0±0.02±0.08 


-0.03±0.04±0.03 


IC 2391 








-0.01±0.02 


-0.12±0.14 


0.01 ±0.02 


-0.07±0.16 


-0.10+0.18 


0.14±0.13 


0±0.02 












-0.04 ± 0.05* 




0.02±0.01 fc 


0±0.02* 
















-0.02 + 0.03 c 












vb 187 


5339 


4.5 


0.92 


0.14±0.04±0.06 


0+0.03 ±0.05 


0.06±0.02±0.08 


0±0.04±0.05 


-0.07±0.04±0.05 


0.12±0.09±0.09 


0.04±0.03±0.05 


iyades Pim i son 








0.13±0.05 


0.01±0.09 


0.05±0.05 


0.07±0.07 


0.03±0.05 






Hyades Fri ei 








0.13 


0.01 


0.04 


0.06 


-0.06 







a Average values were computed by discarding the coolest star R15. b VXR67 and VXR76a were not included in the average. c Average values by considering only stars with T eff > 5500 K 
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Fig. 3. [X/Fe] as function of effective temperatures (T e ff) for both clusters: squares and triangles are IC 2602 and IC 2391 stars, 
respectively. 



field stars by Bodaghee et al. (2003) and by Gilli et al. (2006), 
who explained them as due to NLTE effects. Also in our case 
NLTE is most likely the reason for the different abundance ra- 
tios obtained for warm and cool stars. However, we note that, 
while for Ca we obtain the same decreasing trend as in the liter- 
ature studies, we also find that [Ti/Fe] ratios decrease with T e ft, 
at variance with the increasing trend found by both Bodaghee et 
al. and Gilli et al. We propose that the different behaviour is due 
to the much younger age of our sample stars: they have enhanced 
levels of chromospheric activity and are more likely affected by 
NLTE overionization, rather than by over-recombination; as dis- 
cussed by Takeda et al. (2007), the latter would indeed cause 
enhancement toward lower T e g values. 

To confirm that over-ionization is the reason for the low Ti i 
abundances of cool stars, we computed Ti abundances using a 



few Ti ii included in our spectral range; the results are shown in 
Fig. |4] where [Ti n/Fe]-[Ti i/Fe] is plotted as a function of T e ff. 
The figure shows that the three stars cooler than 5000 K charac- 
terized by lower Ti i abundances indeed have larger differences 
between Ti n and Ti i, confirming the over-ionization hypothesis. 

For Na, Mashonkina et al. (2000) computed standard NLTE 
effects for the lines used in this paper, finding that for logg~ 4.5, 
NLTE corrections are small and do not greatly depend on T e ff. 
Again, we suggest that the low [Na/Fe] ratios that we measure 
for the coolest sample stars, might be due to overionization due 
to chromospheric activity. We also note that, while Ti and Ca 
show lower abundances below ~ 5000 K, the trend of decreas- 
ing Na already starts at ~ 5500 K. Independently of the explana- 
tion for the discrepant abundance ratios for cool stars, we com- 
puted the mean [X/Fe] ratios for Na i, Ca i, and Ti i both con- 
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Table 3. Random errors due to uncertainties in stellar parame- 
ters for one of the coolest stars and for the warmest one in our 
sample. 



SHJM2 


T cff =5970K 


logg=4.45 


f=1.18km/s 


A 


AT cff =±60 


A log g=+0.25 


A£=+0.2 




(K) 


dex 


km/s 


[Fe/H] 


0.05/-0.04 


-0.03/0.04 


-0.02/0.04 


[Na/Fe] 


0.03/-0.02 


0.01/-0.01 


0.02/-0.03 


[Si/Fe] 


-0.04/0.03 


0.03/-0.03 


0.03/-0.03 


[Ca/Fe] 


0.03/-0.01 


-0.03/0.03 


-0.02/0.02 


[Tii/Fe] 


0.01/-0.02 


0.01/-0.01 


0.01/-0.03 


[Tin/Fe] 


0.02/-0.03 


0.08/-0.07 


0.02/-0.01 


[Ni/Fe] 


-0.02/0.01 


0.02/-0.01 


0.01/-0.02 


R95 


T eff =5020K 


logg=4.5 


£=0.82 km/s 


A 


AT eff =±60 


A log g=+0.25 


A<f = +0.2 




(K) 


dex 


km/s 


[Fe/H] 


0.01/-0.02 


-0.01/0.01 


-0.05/0.04 


[Na/Fe] 


0.01/-0.03 


-0.04/0.01 


0.05/-0.03 


[Si/Fe] 


-0.03/0.04 


0.04/-0.02 


0.04/-0.03 


[Ca/Fe] 


0.05/-0.06 


-0.05/0.04 


-0.01/0.01 


[Tii/Fe] 


0.07/-0.05 


-0.02/0.02 


0.01/-0.02 


[Tin/Fe] 


0.05/-0.04 


0.09/-0.08 


0.03/-0.04 


[Ni/Fe] 


-0.03/0.02 


0.05/-0.04 


0.01/-0.01 




4500 5000 5500 6000 



T«„(K) 

Fig. 4. [Ti n/Fe]-[Ti i/Fe] as a function of T e fr for both clusters: 
squares and triangles are IC 2602 and IC 2391 stars, respectively. 



sidering all stars, and discarding the three coolest members (see 
Table|4]i. Average values for Na were also computed considering 
only stars warmer than 5500 K. 

Average values of [Fe/H] and [X/Fe] ratios for all elements 
are shown in Table[4] The average metallicity is [Fe/H]=0+0.01 
and [Fe/H] =-0.0 1+0.02 for IC 2602 and IC 2391, respectively; 
the a-elements, as well as sodium and nickel, also confirm a 
homogeneous solar composition, with no evidence of star-to-star 
scatter. Our results for [Fe/H] are in agreement with those of 
Randich et al. (2001) who, as mentioned, also obtained a solar 
metallicity for both clusters. We stress, however, that the higher 
quality data used here allow us to obtain better measurements 



for each star and thus to place tighter constraints on the average 
and the absence of any star-to-star variation. 

5. Discussion 

As mentioned in the Introduction, along with IC 4665, our IC 
clusters represent the only young PMS clusters for which a de- 
tailed and accurate chemical composition analysis is available so 
far. For both them we obtain a solar metallicity and global abun- 
dance pattern, similarly to the results of Shen et al. (2005) for 
IC 4665. In the following we will discuss the implications of our 
results for the three issues discussed in the introduction. 

5.1. Metallicity and debris disks 

Thanks to the Spitzer Space Telescope, the fraction of debris 
disks has now been determined in a number of open clusters 
and associations with ages up to 600 Myr (the Hyades). These 
observations have allowed investigators to derive the evolution 
with age of the fraction of debris disks. As mentioned in the in- 
troduction, this number reaches a peak around 20-30 Myr and 
both the decay and variation of 24 pm excess ratios around 
FGK-type stars appear similar to those measured around A- 
type stars (Siegler et al. 2007). The question remains, how- 
ever, whether this decay might depend on the metal content. 
Siegler et al. (2007) carried out SPITZER 24 pm observations 
of IC 2391 and found that the fraction of debris disks around 
stars of spectral type B5-A9 is below that measured in other sam- 
ples of similar age and even in older clusters such as the Pleiades 
and NGC 2516. On the other hand, the fraction of debris disks 
around FGK (solar-type) stars is normal. IC 2602 also seems to 
have a lower than expected fraction of debris disks around B and 
A stars. This is shown in Figure 13 of Currie et al. (2008J, al- 
though no reference is provided for S PITZER observations of 
IC 2602. With the caveat that the lower than expected number of 
A stars hosting debris disks could be simply due to a statistical 
deviation, Siegler et al. (2007) suggested possible environment 
effects, among which metallicity, as responsible for this trend. 

To our knowledge, only a minority of the open clusters ob- 
served by SPITZER have an accurate metallicity determina- 
tion available (Pleiades, NGC 2516, Hyades). Vice versa, sev- 
eral clusters in the age range ~ 100 - 600 Myr with metallicity 
measurements do not have SPITZER observations. In order to 
reveal the possible role of metallicity in debris disk frequency 
and evolution, a large sample of clusters with both metallicity 
and SPITZER measurements should be constructed; in particu- 
lar, metallicity should be determined for additional PMS clus- 
ters observed by SPITZER. Nevertheless, the comparison of 
IC 2391 with the Pleiades and NGC 2516 seems to suggest that 
the reason for the lower fraction of A-type stars with 24 pm ex- 
cess in IC 2391 is not metallicity. In fact, both NGC 2516 and 
the Pleiades have a fraction of stars showing an excess in agree- 
ment with the expected decay law from younger clusters and 
field stars, and larger than in IC 2391. Also, both clusters have a 
very close to solar metallicity (see Sect. !5.2l below and Terndrup 
et al. 2002). Our reliable confirmation of a solar metallicity for 
IC 2391 (and IC 2602) thus seems to rule out a metallicity effect. 

5.2. Abundances as tracers of a common origin 

We recall that a common origin has been suggested for sev- 
eral near-by OCs and young associations like the Pleiades, 
NGC 1039, and IC 2602, among others. In order to check 
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whether their chemical composition supports the common ori- 
gin scenario or not, in Fig. [6] the abundance pattern of IC 2602 
is compared with those of Pleiades (both the King et al. 
2000 and Wilden et al. 2002 analyses are considered), and of 
NGC 1039 (Schuler et al. 2003). Note that the two studies 
of the Pleiades yield different metallicities; namely, [Fe/H] = 
0.06 + 0.03 (King et al. 2000), while Wilden et al. (2002) found 
log n(Fe)=7.44+0.02. This translates into [Fe/H]=-0.08, -0.05, 
and -0.1, when assuming our solar Fe abundance, their solar 
model Fe abundance (log n(Fe)=7.49, their Table 3), or their 
adopted value for the Sun (log n(Fe)=7.54 -their Table 5). We 
also mention that the study of Wilden is based on a much large 
number of stars than that of King et al. Given the uncertainties on 
the derived abundances, the comparison between Pleiades and 
IC 2602 does not confirm nor exclude the presence of a common 
origin and a homogeneous analysis is needed. On the contrary, 
the abundance pattern of NGC 1039 appears to be significantly 
different from that of IC 2602, well beyond the observational un- 
certainties, not only for iron, but also for the other elements. In 
particular, NGC 1039 shows an iron content slightly over-solar, 
with A[Fe/H]=+0.07 dex (Schuler et al. 2003) with respect to 
IC 2602; on the other hand, silicon, titanium and nickel seem to 
show a sub solar composition. The [Ca/Fe] ratio is instead close 
to the solar, although in this case the uncertainty is much larger, 
with a standard deviation of ~ 0. 1 dex . 

To investigate if such a behaviour could be due to inhomo- 
geneities in the abundance scales, we re-analyzed one member 
(JP1 33) of Schuler et al. (2003) using their published EWs and a) 
their atomic parameters, b) our atomic parameters for the lines in 
common with our list. Also, we analyzed one of our sample stars 
(R66), adopting their atomic parameters for the lines in com- 
mon. The results of the different tests are listed in Table [5] both 
the [Fe/H] values and [X/Fe] ratios that we obtain for the JP133 
are similar to those derived by Schuler and collaborators; also, 
for R66 we retrieve abundances very similar to our original ones. 
These results show that the two abundance scales are consistent 
with each other. In other words, we confirm that the observed 
discrepancy in abundance pattern is not related to the abundance 
analysis, but most likely reflects the different chemical compo- 
sition of these two clusters. In turn, these different abundance 
patterns do not support a common origin for NGC 1039 and 
IC 2602. 

5.3. IC clusters in the disk 

Based on the fact that the average metallicity of solar-type in 
the solar neighborhood is about 0.2 dex below that of the Sun, 
Wielen et al. d 19961 1 suggested that the Sun might have been 
formed in the inner part of the disk and then migrated to its 
present position. The solar metallicity and solar abundance ratios 
for IC 2602, IC 2391, and IC 4665 do not support this hypothe- 
sis. As mentioned, given their young age, these clusters have had 
no time to move from their birthplace, the solar neighborhood; 
thus their abundances trace the chemical composition of the so- 
lar vicinity at the present time. The fact that they share the same 
composition as the Sun (both iron and other elements) provides 
support to the idea that the Sun was formed where it is now and, 
more in general, that the solar neighborhood underwent very lit- 
tle (if any) chemical evolution in the last ~ 4.5 Gyr. 

In Figure we compare the metallicity distribution of thin 
disk stars from Soubiran & Girard (2005) with the average 
[Fe/H] values of IC 2602, IC 2391, and IC 4665, as well as 
with the distribution of metallicities of all clusters closer than 
500 pc with a good quality measurement of [Fe/H]. The figure 



clearly shows that the metallicity distribution of open clusters 
is different from that of field stars; the last extends to signifi- 
cantly lower metallicities and is characterized by a larger dis- 
persion. The cluster distribution peaks at the same solar metal- 
licity of the three PMS clusters; the mean metallicity of all the 
clusters, including the PMS clusters, is [Fe/H] c i ust . = 0.02 + 0.1, 
to be compared with [Fe/H]g e id = -0.18 ± 0.24 for field stars. 
The fraction of field stars with a very close-to-solar metallic- 
ity (-0.05 <[Fe/H]< 0.05) is 14 % only, while almost 62 % 
have [Fe/H] below -0.1. Also, we note that the oldest cluster 
in the sample (NGC 752, ~ 2 Gyr) shares the solar metallicity 
of the IC clusters. In other words, the inclusion of IC 2602 and 
IC 2391 in the cluster sample not only reinforces the conclusion 
that most clusters in the solar neighborhood are characterized by 
a solar metallicity, but also reveals the lack of an age-metallicity 
relationship. 

The difference between the young cluster and field stars 
[Fe/H] values suggests that a large fraction of field stars have 
not formed where they are now in the solar vicinity; rather they 
originated at larger distances from the Sun, in the outer part of 
the Galaxy, where the metallicity is lower because of the pres- 
ence of a negative radial metallicity gradient (e.g., Magrini et al. 
2009 and references therein). These stars then underwent radial 
mixing (see e.g. Prantzos 2008), moving to their present posi- 
tion. 



6. Summary and conclusions 

In this paper, we present a high-resolution spectroscopic study 
of the two young near-by PMS clusters IC 2602 and IC 2391; 
along with IC 4665, these are the only young OCs whose abun- 
dances have been determined so far. We measured elemental 
abundances of Fe i, Na i, Si i, Ca i, Ti i, Ti n, and Ni i for a 
sample of 15 dwarf members of spectral types G-K. The clus- 
ters show a highly homogeneous solar composition, with no hint 
of dispersion among the members, within the very small uncer- 
tainties. 

Our results strongly favour the hypothesis that the Sun was 
born where it is now, rather than having migrated there from the 
inner disk, and that the solar neighborhood has undergone little 
evolution in the last 4.5 Gyr. On the other hand, the comparison 
with the metallicity distribution of field stars in the solar vicinity 
suggests that metal-poor stars in this sample originated from the 
outer parts of the Galactic disk. 

The comparison between the abundance patterns of IC 2602 
and NGC 1039 seems to exclude a common origin for the two 
clusters, as proposed in the literature; iron as well the overall 
abundance pattern are different, well beyond the measurement 
uncertainties. 

Finally, our results seem to exclude metallicity as the reason 
for the low fraction of debris disks observed among A-type stars 
in IC 2391. We stress however that, in order to derive definitive 
conclusions on the relationship between metallicity and the evo- 
lution of debris disks, accurate metallicity determinations should 
be performed for PMS clusters with available S PITZER obser- 
vations. 
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Table 5. Re-analysis of a NGC1039 member (JP133) and of R66 (IC 2602). 



Star 


T c ff 


logg 


f 


[Fe/H] 


[Si/Fe] 


[Ca/Fe] 


[Ti/Fe] 


[Ni/Fe] 




(K) 




(km/s) 












JP133schulcr 


5710 


4.54 


1.88 


0.07+0.05 


-0.02±0.08 


0.02+0.08 


-0.17±0.08 


-0.13±0.05 


JP133" 


5710 


4.54 


1.88 


0.08±0.06 


-0.09±0.07 


-0.02±0.1 


-0.17+0.11 


-0.16±0.08 


m33 b 


5710 


4.54 


1.88 


0.05±0.14 




-0.02±0.16 




-0.18+0.11 


R66 


5640 


4.45 


1.30 


-0.02±0.07 


-0.02±0.07 


0.05±0.06 


-0.03±0.05 


-0.04±0.05 



11 With the line list by Schuler et al. (2003) and our method 
b EWs by Schuler et al. and our log gf for common lines. 
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Fig. 5. Distribution of [Fe/H] of thin disk stars from Soubiran & Girard (2005). We include only stars in Soubiran & Girard with 
available kinematics and consider as thin disk population stars with P t h.disk ^ 90%. This makes a sample of 358 stars. The average 
metallicities of IC 2602, IC 2391, and IC 4665 are also shown, along with the distribution of [Fe/H] of open clusters within 500 pc 
from the Sun with available good measurement of the metallicity. References for cluster data are: Ford et al. (2005, Blanco 1), Friel 
& Boesgaard (1992, Coma Ber), Paulson et al. (2003, Hyades), Wilden et al. (2002, Pleiades), Boesgaard & Friel (1990, a Persei), 
Sestito et al. (2003, NGC6475), Pace et al. (2008, Praesepe), Schuler et al. (2003, NGC1039), Terndrup et al. (2002, NGC2516), 
Hiinsch et al. (2004, NGC2451a), Jeffries et al. (2002, NGC6633), Jacobson et al. (2007, IC4756), Sestito et al. (2004, NGC752). 
We have considered as thin disk stars those with a probability of membership to the thin disk > 90 %. 
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Fig. 6. Comparison of [X/Fe] ratios as function of [Fe/H] for IC 2602 (our study -star symbol), NGC 1039 (Schuler et al. 2003 
-filled square) and the Pleiades (Wilden et al. 2002; King et al. 2000 -open and filled triangles.) The solid lines indicate the average 
of IC 2602, while dashed lines indicate average [Fe/H] or [X/Fe] +lcr. 



